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快速 射电 暴 观测 数据 干扰 缓解 方法 研究 ” 


刘 艳 玲 1231 陈 卯 蒸 123 
(1 中 


(3 新 疆 微波 技术 相 
摘要 从 海量 的 天 文 观测 数据 中 快速 搜寻 罕见 的 快速 射电 暴 (Fast Radio Burst, FRB) 事 件 , 干扰 缓解 是 其 中 


R (Radio Frequency Interference, RFT) 会 淹没 真实 的 天 文 事件 , 还 会 导致 
其 种 类 的 复杂 性 , 目前 并 没有 一 种 通用 的 方法 可 以 解决 这 


一 项 关键 而 具有 挑战 的 工作 . 射频 干 


搜寻 管线 输出 大 量 的 假 阳 性 候选 体 . 


个 问题 . 为 了 降低 干扰 对 FRB 观 测 搜寻 的 影 
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况 , 针对 主要 的 窄带 干扰 和 宽带 干扰 建立 了 3 
Uh. 将 该 流程 嵌入 到 FRB 色 散 动态 谱 搜寻 (Di 


于 干扰 来 源 及 
H, 分析 和 
恨 次 的 干扰 缓解 处 理 流 程 , 从 而 有 效 缓解 了 观测 数据 的 干扰 污染 情 
spersed Dynamic Spectra Search, DDSS) 管 线 中 , 实验 结果 表明 ， 
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究 了 南山 26 m 射 电 望 远 镜 工 波段 观测 数据 中 的 干扰 情 


搜寻 管线 的 检测 率 和 检测 精度 得 到 了 进一步 的 提高 . 该 方法 为 FRB 观 测 数 据 干扰 缓解 处 理 提供 了 有 价值 的 参考 . 


关键 词 射电 连续 谱 : 暂 现 源 , 射频 干扰 , 方法 : 数据 分 析 
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1 引言 


射电 天 文学 利用 射电 天 文 望远镜 被 动 地 接收 
宇宙 中 极其 微弱 的 射电 信号, 而 接收 到 的 电磁 波 


信号 是 非常 微弱 的 , 大 概 是 150-220 dB- W- m-?2， 


例如 , 最 强 的 射电 源 在 300 MHz 的 功率 谱 密度 大 概 
是 105 Jy (1 Jy = 10-7 W.m-?.Hz-!)H. 因 
射电 天 文 观测 活动 非常 容易 受到 人 为 活动 的 无 线 


此 ， 


干扰, 如 移动 电话 、 了 卫星、 飞行 器 、 广 播 电视 以 


T 


及 台 址 内 电子 设备 等 外. 另 一 方面 , 随 着 射 
技术 的 发 展 , 射电 望远镜 接收 机 系统 灵敏 度 不 断 提 
高 , 带宽 逐渐 增加 , 甚至 超过 射电 天 文 频谱 分 配 保 
PUGH, 而 与 此 同时 , 数字 接收 机 、 数 字 终 端 、 商 
设备 、 电 气 设备 、 ! 线 缆 的 引入 以 及 各 种 航 
空 航天 卫星 通信 业务 的 发 展 使 射电 天 文 台 志 
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磁 环 境 变 得 更 为 复杂 . 在 现在 以 及 未 来 的 射电 天 
文学 中 , 射频 干扰 (Radio Frequency Interference, 
REFI) 绥 解 问题 将 是 主要 的 挑战 之 一 . 单 天 线 观测 最 
易 受 到 干扰 的 影响 , 因为 它们 对 干扰 几乎 没有 识别 
力 , 会 导致 天 文 信号 和 干扰 信号 一 同 进入 系统 刁 . 
快速 射电 暴 (Fast Radio Burst, FRB) 是 一 种 
明亮 的 、 持 续 时 间 仪 为 数 毫 秒 的 射电 天 文 瞬 变现 
RS, 相 较 于 脉冲 星 短 而 稳定 的 脉冲 周期 特性 , 捕 
捉 罕见 的 FRB 就 不 那么 容易 . 在 无 线 电 瞬 变现 象 和 
脉冲 星 搜寻 中 , 理想 化 的 数据 由 弱 天 体 物 理 信和 号 与 
不 相关 的 高 斯 噪声 全 加 组 成 . 而 干扰 的 出 现 让 情况 
变 得 复杂 . 干扰 会 淹没 FRB 信 号 , 降低 信 噪 比 , 导致 
检测 算法 漏 检 天 文 事件 的 发 生 ; 另 一 方面 , 干扰 会 
触发 搜寻 算法 生成 假 的 候选 体 样本 , 增加 后 期 筛选 
困难 . 因此 , 对 于 FRB 的 搜寻 观测 , 干扰 缓解 是 一 件 
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非常 必要 的 工作 . 

射电 天 文 的 观测 研究 一 直 伴随 着 RFI 绥 解 的 难 
题 . 尽管 天 文 台 站 一 般 都 建立 在 远离 城市 的 地 方 ， 
但 是 , 随 着 城市 化 建设 进程 的 推进 以 及 人 类 活动 内 
容 和 范围 的 扩大 , 台 站 的 电磁 环境 只 会 越 来 越 复 
AR. 图 1 展示 了 电磁 干扰 来 源 及 缓解 策略 框架 示意 
图 . 随 着 空间 无 线 电 技术 的 发 展 及 台 站 内 电子 设备 
的 不 断 增 多 , 仅 依 靠 设 立 无 线 电 静 区 、 管 控 台 站 内 
电子 设备 等 管理 方面 的 努力 无 法 免除 射频 干扰 对 
天 文 观测 数据 的 影响 . 因此 , 必须 在 接收 机 链 路 和 
数据 后 处 理 阶段 采取 干扰 缓解 技术 措施 . 在 接收 机 
链 路 中 , 通过 在 射频 模拟 前 端 设计 适当 的 滤波 器 滤 


RFI source 


图 1 电磁 干扰 缓解 策略 框架 示意 图 
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除 可 能 使 放大 器 饱和 的 强 干 扰 信 号 fa. 在 数字 域 
链 路 中 通过 采用 实时 在 线 处 理 和 离线 后 处 理 方法 
对 更 弱 的 干扰 信号 进行 消减 , 具体 措施 包括 时 域 切 
除 B1、 频 域 切 除 岂 和 空间 域 切除 12 于、 自 适应 
WEBS ALAS SUNSETS. 另外 , 在 对 数据 
进行 精细 分 析 时 , 目测 手动 标记 RFI 的 离线 后 处 理 
方式 , 依旧 是 目前 常用 的 方法 . 但 是 随 着 数据 量 的 
增加 , 手动 干预 的 方法 显得 非常 笨拙, 特别 是 对 于 
天 体 物 理 瞬 变现 象 的 搜寻 观测 数据 , 显得 愈加 的 无 
能 为 力 . 总 体 来 说 , 没有 一 种 通用 的 方法 能 够 消除 
所 有 的 干扰 , 在 射电 天 文 领域 干扰 缓解 是 一 项 复杂 
的 工程 化 任务 . 


Fig.1 Schematic diagram of RFI mitigation strategy framework 


随 着 近 几 年 诸如 快速 射电 暴 之 类 的 瞬 态 信号 
观测 研究 工作 的 大 量 开展 , 针对 单 脉冲 观测 数据 的 
干扰 缓解 方法 研究 也 随 之 而 增多 . 但 是 , 由 于 对 实 
时 处 理 的 需求 , 使 得 很 多 基于 峰 度 、 表 面 拟 合 或 者 
傅 里 叶 变换 等 技术 的 干扰 缓解 算法 因为 计算 的 复 
杂 性 而 失效 . 因此 , 对 于 瞬 态 信号 的 处 理 , 考虑 更 多 
的 是 简单 、 快 速 . Men 等 na 提出 Zero-DM matched 
filter (ZDMF) 处 理 方法 进行 瞬时 宽带 干扰 的 消减 
处 理 , 并 同 固定 窗 带 干扰 消除 步骤 一 同 应 用 到 他 
们 自主 开发 的 快速 射电 暴 BEAR (Burst Emission 


Automatic Roger) 搜 寻 管道 中 . Boylet H T 
一 种 边缘 阔 值 分 割 的 干扰 抑制 算法 , 并 与 时 域 阀 值 
切除 结合 , 降低 了 管道 的 FRB 假 阳性 检测 率 . Scloc- 
co 等 PI 开 发 了 开源 的 RFIm (RFI Mitigation) FA 
于 执行 时 域 阔 值 切除 和 频 域 闪 值 切除 , 达到 消除 
瞬 态 观测 数据 中 的 宽带 干扰 和 罕 带 干扰 的 目的 . 
Morello 等 2 开发 了 四 分 位 距 缓 解 (Inter-Quartile 
Range Mitigation, IQRM) 算 法 对 天 体 物理 瞬 态 现 
象 和 脉冲 星 搜 寻 观 测 数据 中 的 窄带 干扰 进行 处 理 . 
总 体 来 说 , 在 对 快速 射电 暴 观 测 数据 进行 干扰 缓解 
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包 暴 观测 数据 干扰 缓解 方 


法 和 
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时 , 会 综合 分 析 观 测 数据 中 的 干扰 情况 , 采取 多 种 
方法 i 以 达到 提高 数据 质量 的 目的 . 
本 文 针对 中 国 科学 院 新 疆 天 文 台南 山 26 m 射 
有 望远镜 FRB 观 测 终 端 产生 的 观测 数据 , 分 析 干 扰 
特点 , 选择 干扰 算法 . 通过 设计 多 层次 的 干扰 缓解 
处 理 策略 对 观测 数据 进行 预 处 理 , 为 FRB 搜 寻 管线 
提供 高 质量 的 观测 数据 , 从 而 提高 搜寻 管线 的 检 出 
率 与 搜寻 精度 . 在 后 续 章 节 中 ， 干扰 分 


方法 . 同时 , 也 是 区 分 来 自 地 面 的 人 为 干扰 和 来 自 
遥远 天 体 的 射电 信和 号 的 主要 特征 之 一 . 


2.2 ”RFI 信 号 特征 及 对 FRB 的 影响 

干扰 按照 其 特点 有 多 种 分 类 描述 方法 , 可 以 是 
瞬时 脉冲 式 的 或 持续 性 的 , 罕 带 的 或 宽带 的 , 周期 
性 的 或 非 周 期 性 的 , 台 址 外 部 的 或 内 部 的 等 . 而 在 
我 们 实际 观测 的 数据 中 , 窄带 持续 性 干扰 与 瞬时 脉 


THA, 将 详细 介 乡 
ETE 


类 及 特点 , 多 层次 的 干扰 缓解 措施 实验 及 应 用 效果 


2E 
=e 


2 ”信号 特征 
2.1 FRB 信 号 特征 

与 脉冲 星 的 单 脉 冲 所 观察 到 的 情况 类 似 , 当 来 
自 遥 远 太 空 的 快速 射电 暴 , 穿越 星际 介质 (Interga- 
lactic Medium, IGM) 被 地 基 望 远 镜 接收 时 , 射电 脉 
冲 会 表现 出 色散 、 散 射 、 闪 烁 等 现象 , 导致 脉冲 轮 
廊 的 展 宽 和 变形 以 及 强度 的 变化 和 频率 的 漂移 . 其 
中 , 色散 的 影响 最 为 显著 , 使 得 其 在 不 同 频率 的 能 
量 信号 到 达 望 远 镜 的 时 间 有 差别 , 具体 表现 为 高 频 
信号 比 低频 信号 先 到 达 , 即 色散 效应 . 由 色散 效应 
引起 的 到 达 时 间 延 迟 表 示 为 : 


all 


At = 4.15 x 10° x (v3? (1) 


H, 频率 v1, v2 分 别 为 两 个 频率 通道 的 中 心 频 

有 位 为 MHz, DM (Dispersion Measure) 为 色散 

ee 分 , 单位 
-pe, 其 表达 式 如 下 : 


—vi?) x DM, 


Mt) 


d 

DM = f nedl , (2) 
其 中 , nw。 是 电子 数 密度 , qd 为 快速 射电 暴 的 距离 , /为 
视 向 路 径 长 度 . 

色散 量 是 天 体 物理 学 中 一 个 重要 的 测量 值 , 通 
过 测量 天 体 信 号 到 达 时 间 随 频率 的 变化 , 可 以 得 到 
色散 量 . 色散 量 越 大 , 代表 爆发 源 越 远 , 通过 色散 量 
可 以 估算 源 的 距离 . 在 脉冲 星 或 者 快速 射电 暴 的 观 
测 技术 中 , 消 色 散 是 最 重要 、 最 具有 特性 的 技术 与 
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冲 式 干扰 又 加 的 现象 是 更 常见 的 . FRB 一 般 为 瞬时 
的 一 次 性 爆发 , 不 可 预知 、 鲜 少 重 复 . 在 目前 观测 
到 的 700 余 例 快速 射电 暴 中 , 重复 射电 暴 只 有 60 多 
Bil, 大 部 分 的 重复 射电 暴 也 只 探测 到 2 次 . 虽然 FRB 
爆发 辐射 的 能 量 很 高 , 但 与 人 为 产生 的 干扰 , 如 天 
空 的 飞机 雷达 信号 、 移 动 通讯 、 卫 星 信 号 及 台 址 1 
子 设备 电磁 辐射 信号 相 比 , 依旧 非常 微弱 , 这 使 得 
FRB 信 号 的 探测 非常 困难 . 目前 针对 FRB 观 测 数 
H, 干扰 缓解 工作 的 难点 主要 在 于 : 1. 干扰 来 源 复 
Ae, 在 观测 数据 上 的 表现 形态 也 是 多 种 多 样 ; 2. F 
扰 是 动态 变化 的 , 不 同时 间 、 不 同 频率 范围 , 干扰 
表现 出 很 大 的 差异 ; 3. 长 期 看 , 干扰 情况 会 持续 恶 
BH. 在 射电 天 文 观测 研究 的 路 上 , 一 直 伴 随 着 干扰 
缓解 方法 的 研究 ,从 未 停止 过 . 针对 FRB 观 测 , 虽 
然 我 们 希望 尽量 减少 数据 中 的 RFI 污 染 情 况 , 但 是 ， 
过 度 切 除 RFI 会 导致 数据 损失 , 甚至 可 能 误 将 天 文 
事件 切除 . 另外 , 在 实时 观测 搜寻 中 , 过 于 复杂 的 干 
爱 解 算 法 或 流程 势必 会 增加 搜寻 管线 的 处 理 速 
此 , 平衡 好 精度 、 速 度 也 是 FRB 搜 寻 观 测 中 
采取 干扰 缓解 措施 所 需要 考虑 的 问题 . 
2.3 ”Peryton 的 特征 

Peryton 即 佩 利 冬 , 传说 中 一 种 半 鹿 半 鸟 的 怪 

S, 而 它 的 影子 却 为 人 形 . 在 射电 天 文 领域 , 用 它 
来 称呼 模仿 快速 射电 暴 的 地 面 无 线 电 干扰 信号 . 
Parkes 64m 射 电 望 远 镜 曾 观测 到 这 类 信号 ， 后 来 被 
鉴定 为 来 自 微波 炉 的 干扰 外 3. 这 类 于 扰 具 有 色散 特 
征 , 很 容易 被 现 有 的 FRB 检 测算 法 误 判 . 但 是 因为 
发 生 的 概率 很 小 , 目前 现 有 的 搜寻 算法 并 没有 对 其 
与 FRB 事 件 做 进一步 的 区 分 . 在 本 文 的 研究 中 , 也 
并 不 对 此 做 进一步 的 研究 和 区 分 . 


已 
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南山 工 波段 信号 干扰 主要 来 源 与 分 
布 


新 疆 天 文 台南 山 基地 始 建 于 1991 年 , 目前 已 拥 
有 26 mm 射电 望远镜 系统 、40 cm 精密 光电 望远镜 系 
统 、1 mm 直径 光学 望远镜 、 多 个 GPS 数据 接收 系统 
和 太阳 色 球 望远镜 等 仪器 设备 , 有 多 台 重 要 天 文 观 
测 设备 正在 建设 或 者 计划 建设 中 . 另外 ,人 类 活动 
的 射电 源 辐射 迅速 扩大 , 可 以 预见 , 未 来 台 址 的 电 
磁 环 境 将 会 更 加 复杂 . 如 图 1 所 示 , 射电 天 文 台 站 
EE 磁 干 扰 来 源 分 为 台 站 内 部 干扰 和 台 站 外 部 干扰 . 
台 站 内 部 干扰 主要 来 自 数 字 时 钟 、 电 源 开关 、 继 
器、 数字 处 理 设备 及 生活 区 电气 设备 (如 无 线 妃 
标 、 电 视 机 顶 盒 、 笔 记 本 电脑 、 微 波 炉 等 ), 其 特点 
是 发 射 信 号 频率 、 时 间 不 固定 , 位 置 相对 固定 . A 
站 外 部 干扰 主要 有 移动 通信 、 电 视 广 播 、 航 空 导 
航 、 航 空 移 动 通信 、 雷 达 、 了 卫星、 业余 无 线 电 等 ， 
其 特点 是 发 射频 段 相 对 固定 , 发 射 功 率 较 大 . 天 文 
观测 数据 中 往往 摊 杂 了 多 种 干扰 源 信息 , 仅 赁 视觉 
上 发 现 的 信号 形式 突变 , 往往 很 难 识别 判断 干扰 来 


3 


a 


识别 
W. 南山 波段 接收 机 主要 用 于 脉冲 星 和 单 脉冲 搜 
寻 观 测 , 因此 , 本 文 主要 针对 工 波段 的 观测 数据 进 
行 FRB 干 扰 缓解 策略 的 分 析 与 研究 . 南山 工 波段 接 
收 机 频率 范围 为 1400-1720 MHz, 在 这 个 频段 范围 
内 的 已 知 外 部 干扰 见 表 1. 


表 1 1400-1720 MHz 频段 范围 内 的 外 部 REFI 来 源 P3l 
Table 1 The sources of external RFI in the 


frequency range of 1400-1720 MHz[23] 


RFI Source Frequency /MHz Characteristic 
GPS 1227.6, 1575.42 Persistent 
Inmarsat 1526-1559 Persistent 
GLONASS 1602 Persistent 
Iridium iiss 9 broad 


bands persistent 


另外 , 从 人 类 社会 活动 时 间 特 点 来 看 , 在 不 同 
的 观测 时 间 段 , 干扰 种 类 、 分 布 、 数 量 以 及 强 弱 也 
是 变化 的 . 新 疆 马 鲁 木 齐 本 地 工作 时 间 是 从 北京 时 
间 上 午 10 点 到 晚上 8 点 , 干扰 状况 处 于 高 度 活跃 的 
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状态 , 从 晚上 8 点 到 第 2 天 上 午 10 点 之 前 , 干扰 相对 


要 少 很 多 . 因此 , 夜间 观测 数据 中 的 干扰 相对 来 说 
会 更 少 一 些 . 


4 干扰 缓解 方法 及 处 理 效 果 
4.1 ”观测 数据 分 析 

针对 FRB 搜 寻 观 测 数 据 进 行 干扰 缓解 , 一 方面 
是 为 了 提高 信 噪 比 , 以 便 让 搜寻 算法 能 够 更 容易 将 
FRB 检 测 出 来 , 另 一 方面 , 希望 通过 减少 干扰 , 降 
低 搜 寻 算 法 误 报 的 概率 . 如 图 2 所 示 , 为 新 疆 天 文 台 
南山 26 mm 射电 望远镜 L 波 段 观 测 数据 的 频率 -时 间 
强度 图 . 由 于 数字 终端 带宽 大 于 工 波段 接收 机 的 工 
作 带 宽 , 因此 , 保存 的 观测 数据 会 有 一 部 分 有 效 观 
测 带 宽 外 信和 号 . 在 实际 的 处 理 过 程 中 , 如 有 需要 可 
以 将 多 余 的 带宽 数据 删除 . 图 中 主要 包含 窄带 持续 
性 干扰 、 多 个 宽带 瞬时 干扰 以 及 短 时 窄带 干扰 . 短 
时 窄带 干扰 在 图 2 中 不 容易 看 出 来 , 但 是 在 本 文 后 
面 的 标准 化 处 理 后 的 降 采 样 数 据 中 会 很 明显 . 本 文 
将 主要 针对 这 3 类 干扰 分 层次 采取 干扰 缓解 措施 . 
干扰 缓解 实验 
窄带 持续 性 干扰 消除 

我 们 选用 的 数据 是 南山 26 m 射 电 望 远 镜 工 波段 
接收 机 观测 的 PSR J1935 + 1616 频 率 - 时 间 二 维 频 
谱 数 据 , 观测 时 间 为 2020 年 5 月 3 日 . 工 波段 接收 机 工 
作 带 宽 为 320 MHz, 数字 终端 频率 带宽 为 512 MHz, 


4.2 
4.2.1 


时 间 分 辨 率 为 64 us, 频率 分 辩 率 为 ] MHz. 在 数据 
处 理 过 程 中 , 仅 保留 了 接收 机 工作 带宽 内 的 数据 . 


将 二 维 频谱 数据 沿 时 间 轴 进行 样本 分 割 , 每 2048 次 
采样 为 一 个 样本 , 得 到 的 样本 为 320 x 2048 像 素 
的 频谱 图 片 . 工 波段 接收 机 工作 带宽 范围 为 1400- 
1720 MHz, 频带 范围 内 出 现 的 窄带 持续 干扰 频段 
范围 有 1450-1465 MHz 和 1535-1545 MHz. XZ Æ 
带 持续 干扰 污染 的 频率 通道 , 通常 的 处 理 方法 是 将 
受 污 染 的 频率 通道 的 数据 直接 删除 . 对 于 删除 的 频 
率 通道 数据 有 多 种 替代 方案 , 如 : (i) 用 固定 常数 进 
行 替代 ; ( 边 用 剩余 通道 或 部 分 通道 中 值 的 均值 蔡 
代 ; GDB RIAT ERRE. 
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图 2” 工 波段 观测 数据 动态 谱 图 
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Fig.2 Dynamic spectrum of L-band observation data 


消 干 扰 在 这 里 只 是 作为 FRB 搜 寻 管 线 的 数据 
预 处 理 过程 之 一 , 考虑 到 后 续 的 过 程 中 , 还 需要 对 
数据 进行 诸如 标准 化 或 者 归 一 化 之 类 的 处 理 , 本 
文选 择 应 用 人 工 高 斯 噪声 进行 替换 . 另外 , 为 了 提 
高 信 噪 比 , 使 数据 具有 更 好 的 显示 度 , 我 们 对 数据 
进行 了 16 倍 的 降 采 样 处 理 . 最 终 得 到 的 待 处 理 样本 
是 320 x 128 像 素 的 频谱 图 . 
4.2.2 ”窄带 短 持续 干扰 消除 

相 较 于 4.2.1 节 的 持续 性 窄带 干扰 , 此 干扰 指 的 
是 只 在 部 分 数据 中 出 现 的 短 持续 性 窄带 干扰 . 对 于 
此 类 干扰 , 应 用 Morello 等 2 提出 的 IQRM 算 法 进 
行 处 理 . IJQRM 的 目标 是 在 给 定 的 动态 谱 数据 块 中 ， 
识别 受 窄 带 干扰 影响 的 所 有 频率 通道 . 首先 是 对 每 
个 频率 通道 数据 进行 频谱 统计 , 我 们 采用 的 是 频率 
通道 标准 差 , 然后 再 通过 计算 比较 确定 哪些 通道 可 
能 受到 了 干扰 污染 . 该 算法 是 非 参 数 的 , 并 且 对 序 
列 数据 中 的 趋势 具有 重 棒 性 . 

频率 通道 标准 差 序列 表示 为 (zi) 呈 7 ， Fea AR 
表 第 i 个 频率 通道 的 标准 差 , n 为 频率 通道 数 . IQR- 
M 算 法 有 两 个 可 以 设置 的 参数 : AEE AE. 半 


径 r 表 示 在 评估 给 定数 据点 的 异常 状态 时 , 设 定 的 
比较 范围 . BERERE EEKE, 代表 高 斯 c (比较 
序列 的 标准 差 ) 的 倍数 , to 为 判断 评估 点 为 异常 值 
的 阔 值 . 具体 步骤 包括 : 

首先 , 通过 输入 动态 谱 数 据 的 频率 通道 标准 差 
序列 , 计算 每 个 频率 通道 与 设 定 半径 范围 内 的 邻近 
通道 标准 差 的 差 值 
AF = T; — Tik. (3) 


a 


RAYE A-r, +r] Lk A 0. 对 于 得 到 的 
序列 A* 用 Qi 表示 , KAKU Tukey? H Fe E 
则 来 识别 其 中 的 异常 值 . 异常 值 的 判断 条 件 如 下 : 


(AF — ur) >tx op, (4) 
其 中 , jx 为 序列 的 中 位 数 , Bur = Qi(0.5), ox 为 序 
列 的 标准 差 , 用 标准 化 四 分 位 距 计 算 如 下 : 


_IQR QQ3— Qi 

1.349 1.349 ， (5) 
其 中 , IQR 为 四 分 位 距 ,@s 表 示 上 四 分 位 数 ,Qi 表 
示 下 四 分 位 数 . 


Ok 
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如 果 (4) 式 的 条 件 成 立 , 则 表示 z; 的 异常 程度 大 ”计算 比较 中 , (4) 式 条 件 成 立 的 次 数 要 少 于 x; 在 所 有 
Fair, 其 为 异常 值 . 但 这 里 并 不 表示 zi 应 标记 为 = 轮 中 被 判 为 高 异常 值 的 次 数 . 
高 异常 值 , 而 zi_# 为 可 以 忽略 的 低 异常 值 , 这 只 是 对 于 标记 出 来 的 频率 通道 , 本 文采 用 与 4.2.1 节 
一 个 相对 的 比较 . 因此 , 需要 综合 多 次 不 同 k 值 计 一致 的 高 斯 噪声 数据 替代 方案 . 图 3 为 两 组 处 理 效 
算出 的 A;, 剔除 其 中 的 低 异 常 值 , 最 后 仅 留 下 高 异 果 图 , 上 图 为 原始 观测 数据 成 图 . 从 图 中 可 以 看 出 ， 
常 值 , 并 对 高 异常 值 所 对 应 的 频率 通道 进行 干扰 缓 ”对 窗 带 干扰 处 理 之 后 , 单 脉冲 的 图 像 质量 得 到 了 很 
解 . 判断 z; 为 高 异常 值 的 条 件 是 : x;_4 在 所 有 轮 的 好 的 提升 . 
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图 3 两 例 罕 带 干扰 缓解 处 理 结果 . 上 图 为 原始 频率 -时 间 图 , 下 图 展示 了 消除 窄带 干扰 后 的 效果 . 


Fig.3 Two examples of narrowband RFI suppression results. The top panels are the raw frequency-time images. The bottom 


panels show the results of mitigating narrowband RFI. 


如 果 在 数据 块 内 的 某 个 频率 通道 上 识别 出 有 固定 窄带 干扰 消除 的 步骤 后 ,再 应 用 IQRM 算 法 

短 持续 性 的 窄带 干扰 存在 , 通过 IQRM 算 法 , 它 就 。 “信号 质量 得 到 了 更 显著 的 提升. 

会 将 所 在 频率 通道 的 数据 进行 干扰 缓解 处 理 . 通常 
3 ”宽带 瞬 乓 缓解 

情况 下 , 窄带 持续 性 干扰 也 会 被 准确 识别 到 , 并 进 si 


4.2. 


行 处 理 . 但 是 我 们 做 了 一 个 实验 , 在 原 图 中 的 1450- 宽带 瞬时 干扰 是 最 为 常见 且 最 昂 引 起 FRB 搜 
1465 MHz 范 围 的 窄带 干扰 , 如 果 我 们 先进 行 高 斯 。 寻 软 件 误 判 的 干扰 , 它 与 PRB 事 件 信号 的 区 别 就 在 
噪声 蔡 换 ,然后 再 进行 IQRM 算 法 处 理 , 对 于 部 分 。 “于 它 不 具有 色散 效应 . 早 在 2009 年 , Batough 等 四 


图 片 , 可 以 进一步 提高 单 脉冲 在 图 片上 的 显示 度 . 提出 了 Zero-DM Filter (ZDF) 方 法 进行 此 类 干扰 的 
图 4 的 上 图 是 未 经 干扰 缓解 处 理 的 原 图 , 中 图 是 只 消除 . 该 方法 通过 对 每 个 数据 减 去 该 数据 所 在 采样 
应 用 了 IQRM 算 法 去 除 窗 带 干 扰 后 的 效果 , 下 图 是 点 的 所 有 频率 通道 数据 的 均值 从 而 达到 消除 干扰 
先 对 1450-1465 MHz 固定 窗 带 干扰 进行 处 理 后 , 再 的 目的 . 此 方法 应 在 消 色散 处 理 之 前 应 用 , 计算 公 
进行 IQRM 算 法 处 理 后 的 效果 图 . 很 显然 , 加 入 对 式 如 下 : 
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7 Penans 致 过 度 切除 信号 的 风险 . Men 等 中 提 出 ZDMF 处 理 

S'i, j) =S S(i, 6 2 ates sas 
= 869) Ta <> CD: (© 方法 , 该 方 法 是 ZDF 的 升级 版 本 , 它 通过 计算 出 每 


其 中 , 5S(i,j) 为 在 第 i 个 频率 通道 , 第 7 个 采样 点 的 数 
HE, ncnans 表 示 频 率 通 道 数 . 5'(i,j) 为 消除 干扰 后 ， 
第 i 个 频率 通道 , 第 ;个 采样 点 的 数据 . 
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图 4 加 入 固定 窄带 干扰 缓解 处 理 的 效果 . 上 图 是 原始 频率 -时 间 图 , 中 
如是 仅 用 IQRM 算 法 去 除 窄带 干扰 后 的 效果 , 下 图 是 加 入 固定 窗 带 干 


扰 消 除 步骤 的 效果 . 


Fig.4 The process effect of adding persistent narrowband 
RFI mitigation. The top panel is the raw frequency-time 
image. The middle panel displays the result of mitigating the 
narrowband RFI with IQRM algrithm only. The bottom 
panel shows the result of adding the step of removing 


persistent narrowband RFI. 


Ht 


里 论 上 , 这 种 方法 可 以 完美 消除 瞬时 宽带 干 
Pi. 但 是 , 在 实际 的 应 用 中 , 会 发 现 宽带 干扰 在 整个 
频率 带宽 上 的 强度 并 不 是 均匀 一 致 的 , 极端 情况 是 
瞬时 宽带 干扰 只 出 现在 部 分 观测 带宽 . 因此 , 会 导 
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个 频率 通道 对 干扰 页 献 的 比例 因子 , 从 而 只 去 除 其 
中 对 干扰 有 贡献 的 部 分 . 图 5 为 分 别 用 两 种 方法 对 
同一 组 数据 的 处 理 效 果 , 其 中 上 图 是 未 经 干扰 处 理 


WER, 中 图 和 下 图 分 别 是 运用 ZDF 和 ZDMF 算 法 
进行 干扰 缓解 处 理 的 结果 , 可 以 看 出 ZDMF 的 效果 
更 好 . 
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图 5 分 别 用 ZDF 和 ZDMF 技 术 对 同一 组 数据 的 干扰 缓解 处 理 效果 . 
上 图 是 原始 频率 -时 间 图 , 中 图 是 用 ZDF 进 行 宽带 干扰 缓解 的 效果 ,下 
是 用 ZDMF 进 行 宽带 干扰 绥 解 的 效果 


Fig.5 The results a 
ZDMF technique respectively. The top panel is the raw 


ter RFI mitigation using ZDF and 


frequency-time image. The middle panel displays the result 
of mitigating broadband RFI with ZDF technique. The 
bottom panel shows the result of mitigating broadband RFI 
with ZDMF technique. 


ZDMF 方 法 的 总 体 思 想 是 通过 将 Zero-DM 波 
形 序列 拟 合 到 给 定 频率 通道 ?序列 的 残 差 最 小 , 得 
到 第 i 个 频率 通道 的 基线 B; 和 比例 因子 ai, 残 差 定义 
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为 : 
x? = (Si — aidamo — Bi), (7) 


其 中 , 5; 是 第 i 个 频率 通道 的 时 间 序 列 . Sg4w_o 是 所 
有 频率 通道 时 间 序 列 的 平均 值 , 即 Zero-DM 波 形 序 
列 . 比例 因子 a 计算 如 下 : 


Sam=0 ` Si 一方 Si È Sam=0 


Qi = (8) 
1 ? 
Sam=0 Sam=o ~~ WN > Sam=0 > Sam=0 
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图 


图 


6 ”两 例 瞬时 宽带 干扰 缓解 结果 . 上 


为 原始 频率 -时 间 图 ， 
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这 里 “表示 为 时 间 序 列 的 内 积 , 和 表示 时 间 序 列 的 


数据 点 数量 . 从 数据 中 去 除 干扰 贡献 部 分 的 计算 公 
式 如 下 : 


S; = Si 一 QiSqm=0 ) (9) 


3; 即 为 干扰 缓解 后 的 第 ;个 频率 通道 的 数据 . 在 处 
理 过 程 中 , 我 们 没有 去 除 偏 置 6;, 因为 它 并 不 影响 


FRB 的 动态 谱 图 样本 检测 . 针对 我 们 的 观测 数据 ， 
它 的 应 用 效果 如 图 6 所 示 . 
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图 , 下 图 是 消除 瞬时 宽带 干扰 后 的 效果 


Fig.6 Two examples of instantaneous broadband RFI mitigation results. The top panels are the raw frequency-time images. 


The bottom panels show the results of mitigating instantaneous broadband RFI. 


最 后 , 为 了 进一步 验证 我 们 设计 采取 3 层 干扰 
缓解 策略 的 必要 性 , 我 们 分 别 单独 使 用 IQRM 
和 ZDMEF 两 种 方法 与 组 合 使 用 IQRM-+ZDMEF 干 扰 
缓解 流程 做 了 补充 对 比 实验 , 实验 效果 如 图 7 和 
图 8 所 示 , 其 中 上 图 是 原始 数据 的 频率 -时 间 动 态 谱 
图 , 中 图 是 仅 采用 IQRM 或 者 ZDMF 处 理 后 的 效果 
图 , 下 图 是 采用 IQRM-+ZDMF 组 合 方式 的 处 理 效 
RE. 可 以 看 出 IQRM-+ZDMF 组 合 的 效果 要 更 优 
于 单独 使 用 IQRM 或 ZDMF 处 理 的 方式 . 


4.3 ”处 理 速度 及 检测 实验 

我 们 将 该 干扰 缓解 处 理 流程 应 用 到 开发 的 
FRB DDSS 处 理 管线 9 中 . 对 一 组 观测 时 长 为 
268.44 s 的 PSR J1935 十 1616 观 测 数据 进行 单 脉冲 
搜寻 处 理 , 运行 结果 如 表 2 所 示 , 加 入 ZDMEF 的 干扰 
缓解 算法 的 检测 流程 所 需要 的 时 间 是 不 加 入 该 组 
解 方法 的 10 倍 以 上 . 这 里 需要 注意 的 是 , 不 同 的 硬 
件 平台 配置 、 程 序 语言 以 及 算法 设计 的 执行 速度 
是 有 很 大 差别 的 . 本 文 提供 的 测试 结果 , 只 是 基于 
目前 的 一 个 单线 程 的 流水 处 理 管线 的 处 理 速度 , 更 
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有 意义 的 是 以 上 两 种 处 理 流程 的 对 比 数 据 . 另外 ， 单 脉冲 数量 , 但 是 这 里 我 们 不 能 用 此 数据 来 进一步 
干扰 处 理 流程 仅 选择 对 窄带 干扰 进行 处 理 时 , 整 计算 和 衡量 搜寻 管线 的 检 出 率 , 因为 脉冲 星 的 每 个 
个 FRB 搜 寻 管线 的 处 理 速度 是 观测 数据 产生 速度 a 甚至 很 多 是 被 淹没 在 观测 
的 6.8 倍 以 上 . 但 就 目前 的 搜寻 管线 而 言 , 如 果 想 要 背景 噪声 中 的 . 


执行 3 层次 的 干扰 缓解 处 理 流程 , 达到 实时 搜寻 的 
目的 , 那么 在 干扰 缓解 处 理 的 环节 上 就 需要 进行 
针对 性 的 设计 , 如 采取 多 线程 的 并 行 处 理 方式 等 
以 满足 速度 上 的 需求 . 从 单 脉冲 检测 数量 上 看 , 加 
入 ZDMF 干 扰 绥 解 算法 的 搜寻 管线 的 单 脉 冲 检 出 


表 2 干扰 缓解 算法 执行 速度 及 FRB DDSS (Dispersed 


Dynamic Spectra Search) 搜 寻 管 线 的 单 脉冲 检测 数 
Table 2 The execution speed of RFI mitigation 


algorithm and pulse detection number produced 
by the FRB DDSS pipeline 


数量 要 更 多 一 些 . 可 见 , 执行 3 层次 的 干扰 缓解 处 
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7 两 例 分 别 用 IQRM 和 IQRM 二 ZDMEF 组 合 的 干扰 缓解 结果 . 上 图 为 原始 频率 -时 间 图 , 中 图 是 用 IQRM 进 行 干扰 缓解 的 效果 图 , 下 图 是 


IQRM+ZDMF44, 


合 进行 干扰 缓解 的 效果 图 


Fig.7 Two examples of the RFI mitigation results using IQRM and IQRM+ZDMF respectively. The top panels are the raw 


frequency-time images. The middle panels display the results of mitigating RFI with IQRM technique. The bottom panels show 
the results of mitigating RFI with combining IQRM and ZDMF techniques. 
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图 8 两 例 分 别 用 ZDMEF 和 IQRM 二 ZDMEF 组 合 的 干扰 缓解 结果 . 上 图 为 原始 频率 -时 间 图 , 中 图 是 用 ZDMF 进 行 干扰 缓解 的 效果 图 , 下 图 是 


IQRM+ZDMF4 FE 


行 干扰 缓解 的 效果 


Fig.8 Two examples of the RFI mitigation results using ZDMF and IQRM + ZDMF respectively. The top panels are the raw 


frequency-time images. The middle panels d 


isplay the results of mitigating RFI with ZDMF technique. The bottom panels show 


the results of mitigating RFI with combining IQRM and ZDMF techniques. 


5 ”总 结 和 展望 
本 文 探讨 了 射电 天 文 观测 所 面临 的 干扰 环境 
挑战 以 及 主要 的 干扰 缓解 措施 , 分 析 了 FRB 信 号 特 


EZ 
IA 


征 及 RFI 对 FRB 搜 寻 观 测 的 影响 , 对 南山 26 mi 
望远镜 台 址 电磁 干扰 环境 、 干 扰 主 要 来 源 以 及 分 
布 做 了 综述 . 为 了 提高 FRB 检 测算 法 的 检 出 精度 ， 
本 文 针 对 南山 LL 波段 观测 数据 中 的 主要 干扰 进行 了 
分 类 , 建立 了 多 层次 的 干扰 绥 解 处 理 流程 , Soe 
RER, 该 流程 取得 了 显著 的 干扰 缓解 效果 . 但 

来 的 另 一 个 问题 是 增加 了 数据 处 理 时 间 , 在 将 该 干 
扰 缓解 流程 直接 应 用 到 FRB 实 时 搜寻 管线 中 时 , 需 
要 提前 综合 考虑 设计 实时 处 理 需 求 . 另外 , 该 方法 
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具有 通用 性 , 可 以 为 其 他 射电 望远镜 的 数据 处 理 


供 有 价值 的 参考 . 
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Radio Frequency Interference Mitigation Methods for Fast Radio 
Burst Observational Data 


LIU Yan-ling'?> CHEN Mao-zheng!?3 YUAN Jian-ping'? LI Jian’? YAN Hao??? 
ZHAI Nan-nan! 


(1 Xinjiang Astronomical Observatory, Chinese Academy of Sciences, Urumqi 830011) 
(2 University of Chinese Academy of Sciences, Beijing 100049) 
(8 Xinjiang Key Laboratory of Microwave Technology, Urumqi 830011) 


Asstract To quickly search for rare fast radio bursts (FRBs) from massive astronomical observation 
data, Radio Frequency Interference (RFI) mitigation is one of the crutial and challenging tasks. RFI will 
produce a large number of false positive candidates, and even submerge real astronomical events. Due 
to the complexity of the sources and types of RFI, there is currently no universal method to solve this 
problem. In order to reduce the impact of RFI on FRB search, the RFI in L band observational data 
of the Nanshan 26 m radio telescope was analyzed and studied. A three-level RFI mitigation procedure 
was established for the main narrowband RFI and broadband RFI, which effectively alleviated the RFI 
pollution of observational data. Embedding this procedure into the FRB DDSS (Dispersion Dynamic 
Spectra Search) pipeline, experimental results show that the detection rate and accuracy of the search 
algorithm have been further improved. This method provides valuable reference for RFI mitigation of 
FRB observational data. 


Key words radio continuum: transients, radio frequency interference (RFI), methods: data analysis 
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